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Q ■ ABSTRACT 

'-^ , We use the near-infrared Br7 hydrogen recombination hne as a reference star 

I ! formation rate (SFR) indicator to test the vahdity and estabhsh the cahbration 

H ' of the Herschel PACS 70 /im emission as a SFR tracer for sub-galactic regions 

52 ! in external galaxies. Br7 offers the double advantage of directly tracing ionizing 

photons and of being relatively insensitive to the effects of dust attenuation. For 

our first experiment, we use archival CFHT Br7 and Ks images of two nearby 

galaxies: NGC 5055 and NGC 6946, which are also part of the Herschel pro- 

"^ ! gram KINGFISH (Key Insights on Nearby Galaxies: a Far-Infrared Survey with 

Herschel). We use the extinction corrected Br7 emission to derive the SFR(70) 

■<sj- ! calibration for Hll regions in these two galaxies. A comparison of the SFR(70) 

^ ! calibrations at different spatial scales, from 200 pc to the size of the whole galaxy, 

reveals that about 50% of the total 70/im emission is due to dust heated by stellar 
populations that are unrelated to the current star formation. We use a simple 

^ I model to qualitatively relate the increase of the SFR(70) calibration coefficient 

j^ I with decreasing region size to the star formation timescale. We provide a cali- 

bration for an unbiased SFR indicator that combines the observed Ha with the 
70 fim emission, also for use in Hll regions. We briefly analyze the PACS 100 
and 160 /im maps and find that longer wavelengths are not as good SFR indica- 
tors as 70/im, in agreement with previous results. We find that the calibrations 
show about 50% difference between the two galaxies, possibly due to effects of 
inclination. 

Subject headings: galaxies:ISM, ISM:structure, infrared: galaxies, infrared:ISM, 
HllI regions 
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Introduction 



Star formation rates (SFRs) have been measured and characterized across all wave- 
lengths, from the X-ray, through the ultraviolet (UV), optical, infrared (IR), and all the 
way to the radio (for a review, see Kennicutt & Evans 2012 for the past decade, and Ken- 
nicutt 1998 for earlier literature). Thanks to the synergy between Spitzer, GALEX and 
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Herschel, many SFR indicators have been recently defined and/or re-calibrated at UV, opti- 
cal, and infrared wavelengths (e.g. Wu et al. 2005; Alonso-Herrero et al. 2006; Calzetti et al. 
2007, 2010; Rieke et al. 2009; Kennicutt et al. 2007, 2009; Boquien et al. 2010; Lawton et al. 
2010; Li et al. 2010; Hao et al. 2011; Murphy et al. 2011a). However, any new calibration or 
re-calibration needs a 'reference ' SFR indicator against which it can be compared. Different 
authors use diverse reference SFR indicators (extinction-corrected Ha, total infrared, radio, 
combination of Ha and 24 /xm emission, etc.); their choices are usually driven by the avail- 
able datasets. Thus, newly (re) calibrated SFR tracers are not usually comparable, to each 
other or to more established indicators. Furthermore, their accuracy and range of applica- 
bility often cannot be reliably established. For example, to calibrate an IR SFR indicator 
with a reference SFR indicator involving another IR band (e.g. the combination of Ha and 
24 fira), degeneracy may arise at high SFR regime where the dust emission dominates. A 
lack of a robust, wide field, reference SFR renders the absolute calibration and universality 
of the 24 fira term in hybrid tracers like Ha+24 /im and FUV-l-24 /im uncertain, even for 
the well-studied SINGS sample (Leroy et al. 2012). Progress in calibrating SFR indicators 
at a variety of wavelengths requires the use of 'unbiased' reference SFR tracers, i.e. those 
insensitive to dust extinction, star formation history, etc. 

Near-infrared (NIR) hydrogen recombination lines, such as Pa (1.8756 /im, from HST, 
Hubble Space Telescope) and Br7 (2.166 /xm), and the radio free-free emission (e.g., with 
GET, the Green Bank Telescope; Murphy et al. 2011a, 2012), trace the ionizing photons 
produced by massive stars and therefore the most recent (<10 Myr) star formation in a 
galaxy. As such, these tracers are less affected by variations in the star formation history 
of a region or galaxy or, as in the case of infrared SFR tracers, by variations in the dust 
content (Kennicutt & Evans 2012; Calzetti 2012). They are also significantly less affected by 
dust extinction than UV and Ha emission; e.g. dust extinction is significantly lower at Br7: 
A_f/Q,=l mag corresponds to ABr7=0.14 mag. The NIR hydrogen recombination lines are 
also free of in-band contamination from adjacent emission lines (e.g., [Nil] for Ha or non- 
thermal emission for the radio free- free emission). Pa and radio free-free emission have both 
previously been used as a reference SFR to calibrate or study other SFR indicators (Calzetti 
et al. 2007; Murphy et al. 2011a). However, the use of these two indicators is limited due 
to the small field of view (FOV) of the Pa observations with HST (e.g. Calzetti et al. 2007 
only studied the central 1-2 kpc of star-forming galaxies), and the time-consuming nature 
of free-free observation to map whole galaxies. 

The WIRCam (Wide Infrared Camera) on CFHT (Canada France Hawaii Telescope) 
provides the capability to observe the NIR hydrogen recombination line, Brackett-7 (Br7, 
2.166 /im, Jones et al. 2002), from the ground with large enough FOV (20'x20') to cover 
most nearby galaxies out to R25 with one pointing (R25 is defined as the 25 mag arcsec"^ 
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isophote). Much like Pa, Br7 emission traces the ionizing photons from massive stars; it 
has the advantage of mapping efficiency over HSTPa or P/3 (Wide Field Camera 3, WFC3) 
observations thanks to the large FOV of modern IR detectors. The lower angular resolution 
provided by ground-based observations relative to the HST is, however, more than sufficient 
for our planned calibration of far-infrared emission. 

High-sensitivity, high-angular resolution infrared telescopes, like the Spitzer and, more 
recently, the Herschel Space Telescopes, have provided a major incentive to investigate dust- 
obscured star formation in galaxies at all cosmic times, and at all scales, from galaxy- 
integrated to sub-galactic scales. Both telescopes have shown that dust-obscured starbursts 
dominate the star formation in the redshift range 1-3 (Le Floc'h et al. 2005; Magnelli et al. 
2009; Elbaz et al. 2011; Murphy et al. 2011b; Reddy et al. 2012). Concurrently, studies of 
nearby galaxies at a few hundred parsec (pc) resolution with Spitzer and G ALEX have yielded 
new insights on the fundamental processes responsible for converting gas into stars and the 
Schmidt Law. Data finally probe the physical scales of Giant Molecular Associations (~100- 
200 pc, Koda et al. 2009), where models of star formation can be discriminated (Kennicutt 
et al. 2007; Leroy et al. 2008; Liu et al. 2011; Rahman et al. 2012; Calzetti et al. 2012). 
Calibrations of SFR indicators in the far-infrared (FIR, >40-50 /xm) regime at sub-galactic 
scales have thus become useful for these kinds of studies, leading recent studies to focus on 
them (Calzetti et al. 2010; Lawton et al. 2010; Li et al. 2010). The integrated FIR emission 
from whole galaxies contains significant diffuse emission that comes from dust heated by 
stellar populations not related to the current star formation. A comparison in Li et al. 
(2010) between the calibrations of SFR(70) for both sub-galactic regions and entire galaxies 
reveals that, on kiloparsec scales and larger, about 40% of emission of 70 /im in the galaxies 
comes from diffuse emission, on average. 

In this paper, we use Br7 maps of two KINGFISH (Key Insights on Nearby Galaxies: a 
Far-Infrared Survey with Herschel, Kennicutt et al. 2011) galaxies, retrieved from the CFHT 
archive: NGC 5055 and NGC 6946, to calibrate the Herschel/PACS 70 /im band as a SFR 
indicator at sub-galactic scales. We also consider the 100 and 160 /xm Herschel bands as SFR 
indicators. The resolution of Herschel at 70 /xm (FWHM~5.5") enables us to resolve regions 
as small as 200 pc in these two galaxies. This scale is considerably smaller than a galaxy, 
and although several Hll regions will be included, we expect the physical characteristics of 
the stellar population in these apertures to be different from those of the entire galaxy. This 
paper is devoted to quantifying the differences of the calibrations between the whole galaxies 
and sub-galactic regions, and to providing interpretation with model comparisons. 

Units in this paper are as follows: ergs • s^^ for luminosity, ergs • s^^ • kpc^^ for lumi- 
nosity surface density (LSD, with symbol S), M© • yr"-"^ for SFR and M© • yr^-*^ • kpc~^ for 



-5- 



star formation rate surface density (S(SFR) or SFRD), unless otherwise specified. The in- 
frared luminosity used in the paper of the IR bands is, L(IR) = z/Lj^, following the common 
definition of monochromatic IR flux. 



2. Data 

Table 1 contains a summary of the data involved in the analyses in this paper, along 
with their angular resolutions. Detailed descriptions are given in following sections. 



2.1. Archive CFHT Data 

Two galaxies, NGC5055 and NGC6946 (Table 2), out of the KINGFISH sample of 61 
galaxies, have archival observations of sufficient depth (>1000 sec integration time in the 
Br7 narrowband, see Section 2.1.3 for the final noise levels of the reduced images) for our 
purposes with WIRCam (Puget et al. 2004), observed by PI Daniel Devost with RUNID 
07AD86 for NGC 5055 and 07BD91 for NGC 6946. WIRCam is the near infrared wide-field 
imager on the CFHT and has been in operation since November 2005. It contains four 
2048x2048 pixel HAWAII2-RG detectors and covers a total 20 arcminx 20 arcmin field-of- 
view with a sampling of 0.3 arcsec per pixel (see Table 2 for the angular sizes of the two 
galaxies) . 

The observations of both galaxies were performed by sequentially shifting the object into 
the center of each detector for each new exposure. Additionally, a certain amount of dithering 
around the object was added to mitigate the effect of bad pixels. This chip-shifting method 
allows adjacent off-target exposures of the same detector to be used for sky subtraction as 
the two galaxies both have angular sizes comparable to the size of one detector. Pipeline 
processed images with sky subtraction (including dark subtraction, fiatfielding with dome 
fiat exposures, basic crosstalk removal, 2MASS absorption estimate and sky subtraction 
using neighboring images) are retrieved from the CFHT science data archive and the total 
number of exposures is listed in Table 3. The actual pixel scale of the retrieved images is 
0.304 arcsec per pixel. Each exposure is stored in a multi-extension fits file with four frames 
corresponding to the four different chips. As each exposure contains only one on-target 
frame, with the off-target frames used for background subtraction, we have a total of 99 
narrowband and 112 broadband sky subtracted frames for NGC 5055 and 61 narrowband 
and 61 broadband sky subtracted frames for NGC 6946. 

In this paper we use "narrowband map" to refer to the observed line-plus-continuum 
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frames obtained with the Br7 filter (central wavelength 2166nni, bandwidth 29.5nm), "Br7 
map" to refer to the Br7 emission line-only frames, and "broadband map" to refer to the 
observed continuum frames obtained with the Ks filter. 



2.1.1. NGC 5055 Reduction 

For the observations of NGC 5055, the different exposures from the same detector are 
well- aligned (as reduced by the CFHT pipeline), but not for exposures from different de- 
tectors. So we first combine all the images from the same detector into one single image 
using IRAF. A sigma clipping algorithm is used to eliminate spurious features, mostly the 
cross-talk from bright stars. All the sky subtracted CFHT images are given a 'chipbias' to 
prevent 16-bit wrapping of negative numbers into high positive numbers and all the bad 
pixels are given a value of 0, so a global 'sky' value for each image, the mode^^ value of the 
whole map, with values close to the 'chipbias', is first subtracted from each image. This 
process reduces our collection of NGC 5055 images from 99/112 (narrowband/broadband) to 
eight: one combined frame for each detector and each filter. The relative astro metry between 
the two bands of the same detector is stable so we use the broadband images to derive the 
alignment transformation between detectors and apply them to both the broadband and 
narrowband images. Thirty point-like sources (foreground stars) are chosen near and within 
the galaxy and the coordinates of their centers are determined. These coordinates are then 
used as input to the IRAF tasks 'geomap' and 'geotran' to align three of the images relative 
to the fourth one, with a manual inspection of the geometrical transformation process to 
mask out outliers due to mis-centering. A higher order fit is used to allow both offsets and 
rotation/stretching to be corrected. The same transformation is applied to align three of 
the narrowband images, relative to the fourth one. The four images from the four different 
detectors are then combined for the narrowband and the broadband separately. 

The chipl and chip4 Ks images of NGC 5055 suffer from severe crosstalk from the 
bright guide star and thus only the chip2 and chip3 Ks images are used for this galaxy. This 
reduced the number of images combined to produce the Ks image to fifty-six 20s exposures. 
However, for the sole purpose of continuum subtraction of the Br7 image, the resultant image 
still provides sufficient depth (see section 2.1.3). 



^^Mode is the global peak of the value distribution. It is robust against the tail of positive sources to 
estimate the sky value. 



-7- 



2.1.2. NGC 6946 Reduction 

For this galaxy we had to proceed with a different approach than for NGC 5055 as the 
different exposures, even from the same detector, is not well-ahgned. No reason is noted on 
the archival images but possibly it is because of the large dither offset from one exposure 
to the other. Thus, for the NGC 6946 images, we manually identified 137 point-like sources 
across the whole field, in order to align all images relative to each other and across the entire 
FOV, also correcting for distortions. The IRAF tasks 'geomap' and 'geotran' are again used, 
with manual inspection, to align the 61 images with reference to the first image for both the 
narrowband and the broadband separately. The rest of the approach remains the same as 
that of NGC 5055 and in this case there is no problem with crosstalk from the guide star. 
Then the 61 images are combined using IRAF task 'imcombine' for the two bands separately. 



2.1.3. Continuum Subtraction and Photometric Calibration 

To obtain the Br7 emission line image, we remove the underlying stellar continuum 
from the narrowband image. We use the broadband Ks image for this purpose. Both the 
narrowband image and the broadband image are first rescaled by the exposure time, and the 
Ks image resolution is matched to the Br7 image resolution with the IRAF task 'psfmatch', 
with a set of point-like sources (~10, these are presumably late type main sequence stars in 
our own Galaxy, and are not saturated in our images). The stars have FWHMs from 1.1" to 
1.4", and the differences between the two bands are less than 0.2". Then, photometry of the 
same set of stars for both NGC 5055 and NGC 6946, respectively, is performed in both bands; 
the average narrowband-to-broadband fractions are determined for both galaxies separately. 
The average fractions serve as a fiducial number we use to perform the continuum subtraction, 
assuming the Br7 features from the stars are negligible. 

Many factors may actually affect the accuracy of this number, such as spectral type 
differences between stars and Hll regions, differential foreground extinction, etc. The Ks 
band targets mostly the Rayleigh- Jeans tail of stellar emission, and, for our galaxies, we 
do not expect large in-band color variations due to variations in stellar populations or due 
to the presence of spatially-variable dust attenuation. Stellar population synthesis models 
(Starburst99) indicate that variations in the stellar flux ratio between the two wavelength 
extremes of the Ks filter are less than 3% for an age range between 6 Myr and 10 Gyr. 
We verify the impact of in-band dust attenuation variations by considering our largest 
Av=3.5 mag derived from H«/Br7. This corresponds to A;^^ ~0.18 mag, when including 
the factor ~2 reduction in attenuation for the stellar continuum relative to the ionized gas 
(Calzetti et al. 2000; Hao et al. 2011). The variations in stellar color across the Ks filter 



due to this value of Aks is about 1%. We also cross-check these scaling numbers against an 
optical continuum subtraction method (Hong et al. 2012), which analyzes the statistics of the 
continuum subtracted images with different continuum fractions and gives a robust estimate 
of the acceptable continuum fraction range. The acceptable range we get from this method 
includes our estimate from rescaling star fluxes, which shows that the two methods agree. 
Having performed this consistency check, we ultimately use the factor determined from 
the stars and inspect the continuum subtracted Br7 images to make sure the subtraction 
looks reasonable. The fractions of the broadband fluxes adopted are 0.0890 and 0.0864 
for NGC5055 and NGC6946 respectively. The contribution of Br7 lines to the Ks filter 
(FWHM 315nm) is negligible (<1.5% even for the youngest star forming regions), thus it 
is not necessary to correct for the over-subtraction of the Br7 line itself when subtracting 
the continuum. In addition to the Br7 line, the H2 rotational transitions may contribute to 
the emission in the Ks filter. These, however, tend to contribute at most 2.4%, and more 
typically less than 1%, to the total Ks band, even in active starburst systems (Calzetti 1997; 
Rosenberg et al. 2013). 

Despite the care taken in the continuum subtraction, we still test for the potential 
effect of incorrect continuum subtraction. For larger the determined value of continuum 
subtraction fraction, we find the smallest multiplicative value that gives over-subtracted Br7 
images for the two galaxies; this value is taken as the upper bound of possible continuum 
subtraction fractions. The symmetric value on the smaller side relative to this is taken 
as the lower bound, as it is more difficult to determine the point where under-subtraction 
occurs. The adopted ranges are 0.0880 to 0.0900 and 0.0852 to 0.0876 for NGC5055 and 
NGC6946 respectively. A set of aperture photometry measurements on Hll regions are 
performed on images with different continuum subtraction levels (the optimal value and the 
two boundaries determined above), which reveals an average of 10%-15% variation in ffux 
measurement between the optimal subtraction fraction and the boundaries. This variation 
depends on how bright the sources are and the difference could rise to nearly a factor of 2 for 
the extremely faint Hll regions, with low signal-to-noise (S/N<1). We thus adopt a S/N>3 
cut that immediately excludes such extreme faint regions (see section 3.3.2). The central 
parts of both galaxies and the centers of the bright stars show artifacts on the line-only Br7 
maps, mainly due to saturation of these regions. The artifacts may also be enhanced by small 
misalignments in the final combined images, because the centroids are poorly determined in 
saturated sources. We make no use of the bright stars and try to avoid the regions (central 
~10") affected by these artifacts in our analysis. Also, because NGC 6946 is in a very crowded 
field, the foreground stars with obvious residual after continuum subtraction are masked 
out. The final Br7 images for NGC 5055 and NGC 6946 both show an uneven background, 
possibly due to an imperfect sky subtraction in the pipeline and the interpolation in the 
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process of aligning the different exposures. To flatten the background, we use the IRAF 
package 'imsurfit' to fit a 2- dimensional low-order surface and subtract the background from 
the images. Although this approach is effective at flattening the line-only images, it may 
not be sufficient for large scale analysis (the peak-to-peak difference in the background for 
the images matched to the resolution of the IR data is about 5 a). We thus avoid deriving 
integrated luminosities for the whole galaxy or radial trends, and limit our analysis to the 
typical (or somewhat larger) sizes of Hll regions. 

The Br7 maps are cahbrated with the header keyword 'phot_cO', which is determined 
using four well modeled spectrophotometric standards (CALSPEC models) and filter trans- 
mission curves, with the information from the WIRCam site^^'^'^'^^. A comparison between 
this calibration and the HST Pa fluxes for these two galaxies, assuming Case-B (Oster- 
brock f989) recombination conditions and similar extinction at Pa (1.8756 fim) and Br7 
(2.166 /im), shows consistency within 5%. The uncertainty of 'phot_cO' is around 0.1 mag- 
nitude (~10%). Including the uncertainty from continuum subtraction, we expect about a 
15% uncertainty in the absolute flux calibration (sum in quadrature of 10% from the cali- 
bration keyword uncertainty and about 10% from the continuum subtraction uncertainty), 
which could be less for bright Hll regions (close to ~10%) as the continuum subtraction 
mainly affects the fainter regions. The flnal continuum subtracted images have l-a noise 
level of 2.12x10"^^ ergs-s~^-cm~^-pixel^^ (0.3 arcsec per pixel) for NGC5055 and 2.48x10"^^ 
ergs-s~^-cm~^-pixel~^ for NGC6946. 



2.2. Herschel PACS Data 

Herschel PACS observations in all three bands (70, 100 and 160 /im) are available for 
both galaxies from the KINGFISH project. The observational program and data processing 
procedures for KINGFISH are described in detail in Kennicutt et al. (2011); Dale et al. 
(2012). Briefly, PACS imaging was obtained with 15'-long scan maps, along two perpendic- 
ular axes for improved image reconstruction, at the medium scan speed of 20"s~^. The 45° 
orientation of the array with respect to the scan direction contributes to a more uniform 
spatial coverage. The raw ('Level 0') data were processed using Version 5 of HIPE (Ott 
2010). Besides the standard pipeline procedures, the conversion from Level to Level 1 data 
included second level deglitching and corrections for any offsets in the detector sub-matrices. 



^■^http://www. cflLt.hawaii.edu/Instrunicnts/Imaging/WIRCam/WIRCamStandardstars. html 
^^http://www.cfht.hawaii.cdu/Instrumcnts/Imaging/WIRCam/dictWIRCam.html=i^UI 



25 



http://www.cilit.hawaii.edu/Instrunicnts/Fihcrs/wircam.html 
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Scanamorphos^^ (Roussel 2012) was used to process the Level 1 PACS scan map data to 
Level 2 final products, where the pixel size is ~ one- fourth of the beam FWHM, i.e., 1.4" at 
70 iim, 1.7" at 100 /xm, and 2.85" at 160 /xm. A multiplicative factor was also implemented to 
update the PACS calibration, which is equivalent to the version 26 calibration^^. We adopt 
the latest values of 10% for all three PACS bands as absolute calibration uncertainty for our 
extended sources^^'^^. The data used in this paper are processed with Scanamorphos version 
12.5 while the most recent KINGFISH data are processed with Scanamorphos version 16.9; 
checks performed against the versions 16.9 of the Herschel maps show <1% difference in 
photometry. 



2.3. Ha Data 

To correct for the small amount of dust extinction in our Br7 maps, we use the Ha 
images of NGC6946 from the Spitzer Infrared Nearby Galaxies Survey (SINGS, Kennicutt 
et al. 2003) and of NGC5055 from the Local Volume Legacy (LVL, Kennicutt et al. 2008; 
Dale et al. 2009; Lee et al. 2011) survey (Fig. 1 &2). The Ha photometry is corrected 
for [Nil] contamination using the [Nil] /Ha ratios listed in Kennicutt et al. (2009, and the 
references therein) (and also in Table 2). The SINGS Ha map of NGC6946 only covers 
the central ~15 kpc by 15 kpc (about 450" by 450") area of the galaxy (Fig. 1), and we 
limit our analysis to this region. We use the Ha map to correct for dust attenuation and 
the typical attenuation corrections of Hll regions are A(Ha)~1.65 mag for NGC5055 and 
A(Ha)~1.12 mag for NGC6946 (see section 2.3.1 and Table 2), after removal of foreground 
Galactic extinction, as derived using the method described in the following section. 



2.3.1. Extinction Correction 

With a simple Case B recombination assumption, the ratio of Ha emission to Br7 
emission from the same Hll region is expected to be ~103 (Osterbrock 1989) for an Hll 
region with electron temperature of 10,000 K and electron density of 10^ cm~^, changing 
by only 35% for electron temperatures in the range 5,000-20,000 K and by only 2%-3% 



^^http://www2. iap.fr/users/roussel/herschel 

^^http://herschel.esac.esa.int/twiki/bin/view/Public/PacsCalTreeHistory?tcmplate=:viewprint 
2*^https://nhscsci.ipac.caltech.cdu/pacs/docs/Photometcr/PICC-NHSC-TR-034.pdf 
^^http://herschel.esac.esa.int/twiki/pub/Public/PacsCalibrationWcb/ExtSrcPhotoni.pdf 
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for electron densities in the range lO^-lO'* cm~^. Figure 5 shows the ratio of S(Ha,obs) to 
E(Br7,obs) as a function of S(70), where the data should lie about the y=103 line if there 
were no extinction. In general, the data for both galaxies show a weak trend for decreasing 
ratio, thus increasing extinction, for increasing S(70). The ratio is mostly below 103, showing 
that extinction is present. Although a couple of regions in NGC6946 have a ratio greater 
than 103, they all have significantly larger error bars indicating lower S/N in at least one of 
the two line measurements. 

We use the deviation of the line ratio from a value of 103 to estimate the extinction of 
both bands. Since A(Ha) = 2.54 • E(B - V) and A(Br7) = 0.36 ■ E(B - V), then we can get 
the nebular gas color excess as 



1 /.poS(Bnsobs)^ 

E(B - V) = - ^ ^(^"'°'^) ' 



0.87 

where E(Br7, obs)/S(HQ;, obs) is the observed surface luminosity ratio. We calculate the 
color excess for each region and apply the extinction correction to S(Br7). The color excess 
is calculated under the assumption of foreground dust extinction (Calzetti et al. 2000), 
because the dust extinction corrections in the NIR regime are relatively insensitive to the dust 
geometry (Calzetti et al. 1996). The extinction corrected E(Br7) is then used as a reference 
SFR indicator. The median E(B-V) values for our regions of NGC 5055 and NGC 6946 at 
70 yum resolution are listed in Table 2, and correspond to a median of 23.6% and 15.7% 
correction of the Br7 luminosity for NGC 5055 and NGC 6946, respectively. 

Hll regions in metal-rich galaxies like NGC 5055 and NGC 6946 often have measured 
electron temperatures that are slightly lower than our default 10,000 K. A temperature of 
8,000 K is more common (Bresolin et al. 2004), which yields an Ha to Br7 ratio of 99.3. For 
this case, the derived E(B-V) from equation (1) would be about 5% lower and the extinction 
correction about 1% lower for the Br7. Given that the conversion from Br7 luminosity 
to SFR(Br7) through Ha luminosity would also be about 4% lower (99.3/103=0.964), the 
SFR(Br7) will be about 5% lower and thus all the derived parameters which are proportional 
to SFR(Br7) in the rest of the paper would be about 5% lower. Since this does not have 
a major impact, we keep our default Ha to Br7 ratio of 103 corresponding to an electron 
temperature of 10,000 K, for generality. 
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3. Analysis 

3.1. Robustness of the Br7 Data 

In order to evaluate our Br7 data reduction and calibration, we derive luminosity func- 
tions for the Hll regions in the two galaxies (section 3.1.1). We also compare a small subset of 
bright Hll regions in NGC6946 with published radio free-free measurements (section 3.1.2). 



3.1.1. Luminosity Functions 

The luminosity functions of the Hll regions in both galaxies are produced from the 
Br7 (uncorrected for extinction) images using the IDL program Hllphot (Thilker et al. 
2000, 2002). Hllphot has an automatic algorithm to identify Hll region seeds and then 
grow the Hll regions from these seeds until it meets another Hll region or reaches a given 
emission measure gradient. Background subtractions are applied to the identified Hll regions, 
and both luminosity and S/N are measured. The luminosity function of the Hll regions 
detected above 5 a are shown for both galaxies, individually and combined, in Figure 3. The 
luminosity functions all have a power law fit, dN/dlogL(Br7) ex; L(Br7)", with power law 
slope a consistent with -1 (shown on Fig. 3), which corresponds to a luminosity function, 
dN/dL(Br7) ocL(Br7)°~^ and agrees with previous studies of Ha luminosity functions (e.g. 
Kennicutt et al. 1989; Thilker et al. 2000, 2002). The most luminous regions have L(Br7) 
~ 10^^'^ ergs-s~^, corresponding to L(Ha) ~ 10^^'^ ergs-s"^ (if no extinction). These values 
are also in agreement with previous findings (Thilker et al. 2000, 2002). This provides a 
consistency check that our continuum subtraction and absolute photometric calibration are 
reasonable. The Hll region luminosity functions have a S-a detection limit of log[L(Br7)] ~ 
35.7, corresponding to star clusters with mass ~ 4 — 7x 10^ M© (STARBURST99, Leitherer et 
al. 1999), assuming an average age of 4-5 Myr and a fully sampled Kroupa (2001) stellar IMF. 
Using the estimates of Cerviho et al. (2002) adapted to our IMF, we expect fiuctuations on 
the ionizing photon rate 6Q{H^)/Q{H^) < 30%. Thus, our regions encompass a sufficiently 
large number of young stars that random fluctuations in the ionizing photon flux are not a 
concern. 



3.1.2. Comparison between Br-f and Free-free Emission 

To further test our Br7 data, we compare the Br7 emission of NGC6946 with the 
thermal free-free component at 33 GHz emission, as observed with the GBT (Murphy et al. 
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2011a, Fig. 4). Murphy et al. (2011a) performed the spectral decomposition of the radio 
and far-infrared emission of the targeted regions, separating the free-free emission from 
other emission components. The GBT observations cover 10 pointings of NGC6946 (Fig. 
1 in Murphy et al. 2011a) and have a FWHM about 25". Our Br7 map of NGC6946 has 
artifacts in the galaxy center, so we only compare the extranuclear measurements. For our 
Br7 emission, we also use an Ha map of NGC 6946 for internal extinction correction (see 
section 2.3 &; 2.3.1), while the radio emission is insensitive to dust attenuation. Both images 
are convolved with a gaussian of FWHM 25" to match the resolution of the GBT 33GHz 
observations. The Enuc 4 (the number four extranuclear region as named by Murphy et al. 
2011a) is the only region among the 10 observed to have statistically significant excess of 
anomalous microwave emission in the 33 GHz band (Murphy et al. 2010); however, this region 
is located outside the coverage of our Ha map, implying that we do not have extinction- 
corrected Br7 flux for this region. Thus, our comparison between the Br7 and the free-free 
emission uses a total of eight regions. A test with the extinction correction for the regions 
before and after the convolution shows less than 5% difference. As our analysis is on a 
region by region basis, we adopt the extinction correction after the convolution for our Br7 
measurements. Photometry is measured with apertures and background annuli similar to 
those used by Murphy et al. (2011a). 

The galaxy-wide star formation rate calibrations do not apply precisely to a single- 
age young populations like those in Hll regions. Here, our apertures correspond to regions 
consisting of a few Hll regions, and both the free-free and the Br7 measurements are actually 
reflections of the ionizing photons, Q(H°), from the young stars. However, for a more direct 
comparison with previous results, we still calculate the SFR(Br7) using the estimated Q(H°) 
flux (Osterbrock 1989), with extinction corrected Br7. SFR"^(33GHz), the star formation 
rate derived from the thermal free-free component of the 33 GHz radio emission, is calculated 
using Equation (11) from Murphy et al. (2011a). Both quantities are calculated assuming 
electron temperatures of 5000K (triangles), lOOOOK (circles) and 20000K (squares) to account 
for the range of variations in electron temperatures found in Hll regions, which depend 
on galacto-centric radius and metallicity (Shaver et al. 1983). The electron temperature 
directly affects the intrinsic Ha/Br7 ratio (Osterbrock 1989), thus affecting the extinction 
correction for the Br7 emission using the Ha emission. It also affects the conversion both 
from measured luminosities of Br7 and 33GHz to Q(H°) flux and from Q(H'') flux to SFR, 
as shown in Equation (11) from Murphy et al. (2011a). Values calculated using Br7 emission 
without internal extinction correction, assuming a typical electron temperature of lOOOOK, 
are shown as open diamonds, representing about 3%-18% corrections for the fllled circles in 
the Br7 emission (with a median ~12%) after removal of foreground Galactic extinction. 
The two inferred SFRs show about 27% difference on average, for the preferred temperature 



-14- 



as shown with the filled symbols. Even when the extreme variations in electron temperature 
(from 5000K to 20000K) are taken into account, the two inferred SFRs are still in good 
agreement. This provides again a consistency check that our data reduction and calibration 
are reliable. Furthermore, our comparison shows that the Br7 data, once corrected for dust 
extinction using the Ha, do not appear to require additional dust corrections, implying that 
the fraction of ionizing photons completely obscured by dust at 2.16 /im is negligible. 

Our 3-0" detection limit within the 12.5" radius aperture is SFR(Br7)~0.0017, or about 
one order of magnitude fainter than the faintest region analyzed by Murphy et al. (2011a). 
The higher sensitivity of Br7 imaging suggests that this may become a preferred approach 
to derive a dust-free SFR indicator over extended sections of star-forming galaxies. 



3.2. PSF-Matching 

In order to measure the photometry consistently among the different maps (Ha, Br7 
and 70 /xm), we first need to convolve both the Ha and the Br7 maps to the resolution of 
the 70 /xm map (~ 5.5"). The point-spread function (PSF) FWHM of the Br7 is about 1.1" 
and the Ha only slightly better (about 1"). 

For the PACS PSFs, we choose the observed PSFs of the asteroid Vesta available on the 
PACS website^°. We use the PSF with the same scan speed (20"s~^) and similar orientation 
(42 deg) to the KINGFISH observations. We measure the curve of growth of the observed 
PSF and find that more than 2/3 of total fiux is included within a radius of 6", slightly larger 
than the PSF FWHMs, 5.5". Although we think this specific observed PSF should be the 
best match to our observations, the curve of growth shows less than 1% difference among 
observed PSFs with different observational parameters. The model PSFs show a different 
curve of growth but they still contain more than 2/3 of total flux within the radius of 6". 
Thus we choose 6" as our aperture radius for the 70 /xm resolution. The PACS PSFs are not 
spherically symmetric and are dependent on the scan direction, thus the PSFs are rotated by 
the scan angle (relative to north, as our maps are reprojected with north-up) of our PACS 
maps. A visual comparison of the rotated PSF and the brightest point-like sources on the 
FIR maps of both galaxies shows differences of only a few degrees of rotation. A comparison 
between the photometry of Br7 maps convolved with one PSF and the same PSF rotated 
by 60 degrees shows only about 7% difference. 

Finally, the Ha and Br7 maps are convolved with the observed PSF, rotated with the 



30 
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scan angle of the Herschel observation for each galaxy, to produce the Ha and Br7 maps at 
70 fira resolution. The convolved Ha and Br7 maps are then registered to the same pixel 
sizes as the 70 /xm maps with astrometry matched using the brightest sources in the galaxies. 
Although the PSF FWHM of the Br7 is small compared to the 70 /xm, we still convolve the 
70 /xm images with the extracted PSF from several unsaturated bright stars in the CFHT 
narrowband images to better match the resolution between the two. We are confident that 
the convolution brings the optical, NIR and 70 fira maps to similar resolution and adds a 
photometric uncertainty no larger than 2%, combining both the difference in observed PSFs 
with different observation parameters and possible small mismatch of PSF rotation angle. 
The final PSF FWHM for the 70 /xm images, the Br7 images and the Ha images is ~5.6". 



3.3. Photometry in Selected Regions 

3.3.1. Aperture Photometry 

Star forming regions are selected on the Br7 images at the original resolution, by placing 
apertures with radii of 6" centered on the peaks of ionized gas emission. Care is taken so 
that apertures do not significantly (<5%) overlap. Further inspection of the 70 /im maps 
and the degraded Br7 and Ha maps is performed to make sure that each aperture covers 
at least one 70 /xm emission peak and is as close as possible to the emission peak both at 
70 fira and in the emission line images. Sources near the continuum subtraction residuals of 
bright stars and the galaxy center are not included as to avoid possible contamination from 
the residuals. The central region of NGC 5055 is identified to have an AGN by Moustakas et 
al. (2010), and it shows one single bright source in the center covering up a region of about 
10" radius (~ 400pc). This region is also avoided when selecting apertures. We selected 
a total of 73 regions for NGC 5055 and 94 regions for NGC 6946 (Fig. 1 & 2). We note 
that these regions are not the same as those described in section 3.1.1; Hllphot could not be 
effectively applied to our resolution degraded images, and we resorted to a manual selection. 
In addition, due to the lower resolution, the regions used for the bulk of our present analysis 
are larger (and more luminous) than those used for deriving the luminosity functions. 

Following the studies of Calzetti et al. (2005), Li et al. (2010) and Liu et al. (2011), 
local background subtraction is essential for photometric measurements in spatially resolved 
regions in order to minimize the effect of diffuse IR emission that is not related to current star 
formation activity. Furthermore, it will serve to mitigate the effects of the uneven background 
in our Br7 maps. We perform the local background subtraction with two methods: one is 
to identify a local background region for different apertures (Calzetti et al. 2007; Li et al. 
2010); the other is to use an annulus outside each aperture not much larger than the aperture 
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radius itself, a method that follows standard practice for photometric measurements. For the 
first method, we choose two local background regions within each of the two galaxies, with 
relatively even background on the Br7 map and similar physical conditions, i.e. on the same 
spiral arm, to ensure that the diffuse IR emission is similar, by visual inspection. For the 
second method, we adopt an annulus whose inner radius and outer radius are two times and 
four times the aperture radius, respectively. The width is chosen to be larger than two times 
the PSF FWHM and is outside of the second Airy ring. Each annulus is thus close enough to 
its aperture to minimize the effects of the uneven background while still containing enough 
pixels for a good local background determination and not including significant emission 
from the region itself. For both methods, the mode is used for the background level due 
to crowding in some regions. The background levels between the two methods agree with 
each other with differences of less than several percent when comparisons are performed in 
uncrowded regions. The results prove that both methods are valid for effectively removing 
the local background in our maps. Due to the uneven background in our Bvy maps, we 
choose the annulus method, which is more robust against such variations. 

The 70 fira resolution aperture correction from a 6" radius to infinity is 1.466, which 
is determined using the PACS observed PSFs convolved with the Br7 PSF. The local back- 
ground subtraction method described above is included in the aperture correction. We use 
aperture corrections to best recover the true luminosity of each source we select, but they do 
not affect the analysis, as all the images are degraded to the same resolution. The aperture 
corrections are only perfectly valid for point sources; however, our sources are usually com- 
pact and the underlying extended emission is relatively weak, thus this correction provides 
a reasonable estimate of the true luminosity of our sources. 



3.3.2. Error Estimate and Signal-to-noise Criteria 

We correct the Ha and Br7 photometry for foreground Galactic extinction (Schlegel et 
al. 1998; O'Donnell 1994) with a single value (Table 2), which is significant for NGC6946 as 
it is close to the Galactic plane, while the Galactic extinction for the FIR bands is negligible. 
This correction could be a potential source of systematic error if the adopted values are not 
correct. 

The major uncertainty in the photometry is the calibration uncertainty: 10% at 70 /im, 
15% at Br7, combining both the absolute calibration uncertainty and the continuum subtrac- 
tion uncertainty, and 10% at Ha (see Section 2). The uncertainty in the local background 
level is also included. We take a 2% uncertainty in the convolution with the observed PSF 
at different rotation angles for the three bands. Other uncertainties, such as the global 
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background uncertainty, are relatively small. The final photometric uncertainty is usually 
dominated by the calibration uncertainty; the uncertainty in local background estimate af- 
fects the faint regions more than brighter regions. 

We accept as 'good data' regions with S/N>3 for all the three images at 70 /im, Ha and 
Br7. The noise is determined using the variation around the mode in the local background an- 
nuli of each region. This S/N cut removes the faint regions whose Br7 photometry is compli- 
cated by the uneven background and uncertainties in the continuum subtraction. Our S/N>3 
cut corresponds to different luminosity cut-offs for different regions, as some of the regions fall 
in high background noise areas. In our 'good data', only five regions have S/N<5, which are 
regions in the East and West outskirts of the galaxies (Figures 2 and 1). About 77% of the re- 
gions have S/N greater than 10. The least luminous region has a L(Br7)=10^^'^ erg-s~^ with 
S/N=10.7, corresponding to a S/N=3 limit of L(Br7) = 10^^'^ erg-s^^ in the low background 
noise area. The lowest S/N region (S/N=3.5) has a L(Br7)=10^^'''' erg-s^"*^, corresponding to 
a S/N=3 limit of L(Br7)=10^^'^ erg-s^^ in the high background noise area. We note that 
the minimum L(Br7)= 10^^-^ erg-s^^ is about eight times higher than the limit discussed 
in section 3.1.1 for the Br7 luminosity function. This difference is explained by the larger 
apertures we use in this analysis relative to those of section 3.1.1, which causes blending of 
Hll regions. The larger apertures are dictated by the larger 70 /im PSF than the native Br7 
PSF. Thus, the same argument about the random sampling of the stellar IMF not being an 
issue when deriving SFRs applies to our overall analysis, as well. Comparing the luminosity 
distribution of our selected regions with the luminosity functions derived using Hllphot (see 
section 3.1.1), we find a 50% completeness limit at L(Br7)~10^^ erg-s~^. 80% of the regions 
we use for the L(70)-SFR analysis have luminosity brighter than this limit, implying that 
our results are robust relative to completeness issues. For the 'good data', we have median 
uncertainties (and uncertainty ranges) of about 17% (15-33%) in Br7, 11% (10-30%) in Ha, 
10% (10-19%) in 70 /xm. The final number of regions used in our analysis is 68 for NGC5055, 
and 92 for NGC 6946. 



3.3.3. Other Sources of Bias and Uncertainty 

Absorption of Lyman continuum photons (LyC) by dust is a potential source of bias 
that decreases the luminosity of hydrogen recombination lines and of the free-free continuum. 
Estimating this effect observationally is challenging, and models are usually used to estimate 
the impact of this effect. We follow the parametrization of Dopita et al. (2003), where the 
lost fraction of H/9 luminosity to LyC absorption by dust is modeled as the product of 
the oxygen abundance and ionization parameter. For the typical oxygen abundance (about 
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solar) and ionization parameter range (U~10^^) expected in our regions, we estimate the 
LyC absorption by dust to be at the level of 10% or less. Hunt & Hirashita (2009); Draine 
(2011) suggests similar scenarios except in very extreme conditions. We thus neglect this 
small negative contribution in our analysis. 

Leakage of ionizing photons out of the regions is also a potential source of bias. Many 
studies (e.g. Ferguson et al. 1996; Thilker et al. 2002; Oey et al. 2007) place the fraction 
of ionizing photons escaping from Hll regions at the 40%-60% level, as evaluated from the 
observed Ha. When accounting for the differential dust attenuation between Hll regions and 
the diffuse medium, the fraction of diffuse Ha is reduced by roughly a factor 2 (Crocker 
et al. 2012). We thus expect that our Br7 measurements will be systematically biased low 
compared to the number of ionizing photons produced by about 20-30%. The dependence of 
this bias on the region's characteristics (luminosity, density, etc.) is still under investigation, 
with different conclusions reached by different authors (Pellegrini et al. 2012; Crocker et al. 
2012). Thus, although the systematic ~20-30% bias is likely to be present for all or most of 
our regions, we will not attempt to correct for it, in light of the uncertainty just discussed. 



4. Results and Discussions 

4.1. SFR(Ha+70) 

Calzetti et al. (2007) and Kennicutt et al. (2007, 2009) proposed to use a combination of 
24 /im emission and Ha emission as an 'unbiased' SFR indicator, where the 24 /xm emission is 
used to represent the attenuated SFR, i.e. attenuated Ha emission. Following these studies, 
we compare the extinction correction and the 70 /xm emission, to see if a similar combination, 

S(Ha, corr) = S(Ha, obs) + aS(70) (2) 

can be established using the 70 /xm emission instead of the 24 /xm emission, where corr stands 
for extinction corrected. Figure 6 shows S(70)/S(Ha,obs) as a function of lo*^-"^^*^^") — 1, 

since 

_ E(Ha,corr) - S(Ha,obs) _ (100-4A(h«) _ i)E(Ha,obs) 

" ~ S(70) ~ S(70) ■ ^^^ 

Linear fits in log-log space, log[S(70)/S(Ha, obs)] = ao + boTog[10°'*^'^*^^"-' — 1] and linear 
fits with fixed unity slope in log-log space, log[S(70)/S(Ha,obs)] = ai + log[10°-''^*^^") — 1], 
are shown as dotted lines and solid lines respectively and the fit parameters are listed in 
Table 4. The best-fit slope of NGC 6946 is consistent with unity and if we adopt the forced 
fit with unity slope, we obtain a = I/IO'^^ which is also listed in Table 4. The l-a dispersions 
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of the data around the hnear fit are also hsted; these dispersions provide a better measure of 
the accuracy of the hnear fit, as the formal uncertainties listed in Table 4 (and subsequent 
Tables) reflect only the internal dispersion resulting from the fit procedure. We use the IDL 
routine 'fitexy' to perform the fit process here and for later fits. This takes into account the 
uncertainties in both axes and derives a straight line fit using ^-square minimization. 

We further test the robustness of our fits, as it seems that the linear fit is determined 
by the highest four data points and thus larger uncertainty may occur without those data 
points, for the NGC6946 fit. Though there is no physical reason to not include these points, 
we obtain a fit slope of 1.062(±0.111) if we exclude them from the fit, which is not too much 
different from the original fit slope 1.026(±0.077) (Table 4) and the fit uncertainty is only 
slightly increased. Thus, it reveals that our fit is actually robust. 

The linear trend obtained for NGC 6946 indicates that using the PACS 70 fxra band as 
a complementary SFR tracer to correct for the attenuated Ha emission is a viable method, 
with a corresponding a value of 0.011(±0.001). We obtain a best-fit slope significantly 
smaller than unity for NGC 5055, probably due to its larger inclination (see Section 4.4). To 
test the apphcability of the derived a, we plot log[S(HQ;, obs) + 0.011S(70)] as a function of 
log[E(HQ;, corr)], the extinction corrected Ha luminosity surface density (Fig. 7), for these 
two galaxies. The general trend of the data is close to the 1-to-l correlation, indicating 
that a = 0.011 works on average. However, the NGC 5055 data lie systematically below 
the 1-to-l correlation. This deviation is discussed in more detail in the next section. A 
comparison of this derived a = 0.011 to the proportionality parameter for Ha + a • 24 /xm 
(a = 0.031, Calzetti et al. 2007), gives a ratio L(70)/L(24)~3, which is consistent with the 
model predictions by Draine & Li (2007) with an average stellar radiation field strength U 
about 300, the value estimated for star-forming regions (Calzetti et al. 2007). The current 
sample size is too small to enable a detailed comparison between the 24 micron and 70 micron 
mixed indicators, which we defer to a later paper where a larger number of galaxies will be 
analyzed. 



4.2. SFR(70) 

In section 4.1 and in this section we seek to obtain a linear trend between the sum of 
Ha plus the 70 /im emission and Br7 and between the 70 fira emission and Br7, respectively. 
At first sight, it will appear that seeking linear trends (slope of unity in a log-log plot) in 
both cases is contradictory. However, we note that the range of luminosities for which each 
calibration is valid will be different and only partially overlapping (see, e.g., discussion in 
Kennicutt et al. 2009; Calzetti et al. 2010). The 70 /im emission will generally be stronger 
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in dust obscured regions where Ha emission is weak. 

We express the IR (70, 100 and 160 /im bands) SFR cahbrations in this paper as 

SFR(IR band) (Mq ■ yr-^) = C^ band.region X 10-^3 X L(IR band) (ergs • s"^) (4) 

We refer to the C as the Cahbration Coefficient. Thus, the Cahbration Coefficients in Calzetti 
et al. (2010) and Li et ah (2010) are Cy^gaiaxy = 0.58, and C7o,7oopc = 0.94 respectively. Using 
the 70 /im to the bolometric infrared luminosity ratio, Lawton et al. (2010) has derived 
C7o,io-300pc = 0.97 for the Magellanic Clouds. 

Adopting SFR = 5.45 x lO'^^ . L^^a) (Calzetti et al. 2010) and L(Ha) = 103 ■ L(Br7), 
we derive SFRs for our regions using their extinction corrected Br7 luminosities. Figure 
8 shows the E(70) as a function of S(SFR). The linear ffis in log-log space, log[E(70)] = 
ao + bo Tog[S(SFR)], and the fixed unity slope fits, log[S(70)] = ai + log[S(SFR)], are shown 
as dotted lines and solid lines respectively and the fit parameters are listed in Table 5. The 
1-0" dispersions of the data around the linear fit are also listed. The fits all have slopes close to 
unity and so we adopt the fit results with unity slope to obtain the Calibration Coefficients, 
C = 10*^~"-^, used in Equation (4). These are also hsted in Table 5; the C7o,2oopc = 1.18(0.02) 
for the combined data (at ~200 pc scale) will be used for the comparison with the Calibration 
Coefficients, C70, from Calzetti et al. (2010) and Li et al. (2010), as larger samples of galaxies 
were combined to derive the calibration in those works. 

The calibration coefficient of NGC 5055 is higher in value than that of NGC 6946 (Table 
4). We provide analysis with the longer wavelengths and reasons in Section 4.4 that lead 
us to believe that the anomalous behavior is in NGC 5055. We suggest that the higher 
inclination of NGC 5055 relative to NGC 6946 prevents an effective removal of the diffuse 
dust along the line of sight to the Hll regions of the galaxy and/or causes over-subtraction 
of the star forming regions contributing to the local 70 micron emission. Although we use 
circular annuli around each region to remove the foreground/background diffuse component 
of dust emission, the high inclination of NGC 5055 implies that there is more intervening dust 
emission within each photometric aperture that is not directly related to the Hll region(s), 
and removal of such a component is subject to larger systematic uncertainties. Indeed, as we 
show in section 4.4, going to lower resolution (longer wavelengths) we increase line-of-sight 
confusion and the NGC 5055 data further deviate from those of NGC 6946. Draine et al. 
(2007) show that the mean interstellar radiation field is 2.5 times weaker in NGC 5055 than 
in NGC 6946, which also manifests itself as a global infrared spectral energy distribution 
(SED) that has a slightly colder mean temperature (Skibba et al. 2011; Dale et al. 2012). 
Thus, the intervening diffuse component in the NGC 5055 apertures is made of dust that is 
slightly colder than that in NGC 6946, which accounts for the behavior observed in section 
4.4 at 100 and 160 fira. In addition to the possible over-subtraction, this also helps account 
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for the higher cahbration coefficient in NGC5055 than in NGC6946: the peak of the diffuse 
dust emission in the former galaxy is shifted to longer wavelengths than in the latter. In fact, 
Hll regions in NGC5055 show higher extinction values than in NGC6946 (Figure 5), which 
may imply higher self-shielding for the dust surrounding the Hll regions, resulting in cooler 
mean dust temperatures (e.g. Calzetti et al. 2010). Both effects are likely to contribute to 
the -50% higher C7o,2oopc in NGC 5055. 

We also revisit the data for both NGC 5055 and NGC 6946 in Li et al. (2010) to un- 
derstand the ~50% higher C7o,2oopc of NGC 5055. We find that the NGC 5055 data in that 
work already show a 2-0" deviation from the mean trend in the sense of being systematically 
higher than the mean, while the NGC 6946 data lie about the mean trend. Indeed, the 
corresponding Cyojoopc value of NGC 5055 is also about 50% larger that that of NGC 6946 in 
the data of Li et al. (2010). This strongly suggests that a larger sample with more galaxies 
is needed to average out fluctuations and derive a more reliable general calibration. 

Although we find roughly linear trends both in section 4.1 and this section, we should 
be reminded that our study only involves two galaxies, and uncertainties are large. A study 
involving a larger sample (which we are planning in the near future) will be better suited 
for analyzing similarities and differences in the two calibrations. 



4.3. SFR Calibrations at Different Scales 

Li et al. (2010) found, from a comparison of the Calibration Coefficients between 70 /im 
emission in sub-galactic regions (Crojoopc) and for integrated (star-forming and starburst) 
galaxies (Cyo^gaiaxy); a ~40% difference between the two values, which they attributed to 
dust heated by stellar populations not related to the current star formation activity. The 
comparison is performed using the following steps: for the calibration in Calzetti et al. (2010), 
a galaxy with a SFR of 1 MQ-yr"^ implies a 70 /im luminosity of 1.725 x 10^^ ergs-s~^; the 
calibration of Li et al. (2010) shows that ~700 pc sub-galactic star forming regions with the 
same total SFR of 1 MQ-yr~^ have a total 70 /xm luminosity of only 1.067 x 10^^ ergs-s~^; thus, 
the difference in these two calibrations reveals an average of ~40% excess 70 /im emission 
in the galaxies on scales larger than 700 pc. With the results in Table 5, we add another 
C70, derived from the two galaxies, with an average region size of ~210 pc (Fig. 9). C70 
increases as we probe star forming regions with smaller size, which gives a factor 2 difference 
between the 200 pc calibration and the whole galaxy calibration, implying that ~50% of the 
70 ^m. emission beyond the 200 pc scale is due to heating of dust by stars not associated 
with current star formation. 
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In the comparison above, we assume that the dust heated by young stellar populations 
in whole galaxies is a scaled-up version of dust heating in 200 pc regions, as traced by the 
ionizing photon flux. This is equivalent to assuming the UV photons leaked from the Hll 
regions trace the leakage of ionizing photons, and that this ratio is constant over all spatial 
scales. The leakage of ionizing photons out of Hll regions is not included in our estimate of 
the S(SFR). If we were to systematically correct all the Br7 luminosities by ~30% (section 
3.3.3), the estimate of the diffuse fraction as inferred from the comparison between C7o,2oopc 
and C-jQ^gaiaxies would increase by a similar percentage, from ~50% to ~65%. This new 
estimate, however, does not include the leakage of dust heating UV photons from the Hll 
regions; if we assume, as done so far, that they trace the leakage of ionizing photons, the 
estimate of the diffuse 70 /im fraction is reset back to ~50%. 

Lawton et al. (2010) provide estimates for the size of Hll regions in the JR. Their analysis 
yields mean 70 /im radii of 60(±20) pc for Hll regions in the LMC and 80(±30) pc for Hll 
regions in the SMC. The radii are calculated as encompassing 95% of the light. Our ~200 
pc radii apertures are thus likely to include most of the UV photons responsible for heating 
the dust within and around Hll regions; the 50% estimate for the diffuse fraction outside 
200 pc radii is thus likely to be a conservative one. 

The interpretation we give for the trend in Figure 9 is that smaller regions contain 
smaller fractions of diffuse IR emission heated by stellar populations that is not related to 
current star formation, i.e. populations older than about 10 Myr. We make use of a simple 
model (Li et al. 2010) to try to qualitatively support this interpretation. The simple model 
involves using, a) STARBURST99 to produce a stellar spectral energy distribution for a 
given star formation history; b) the attenuation curve from Calzetti et al. (2000) applied 
to the stellar light to produce the total IR emission; c) the dust model from Draine & Li 
(2007) to predict IR emission in a given band (here 70 /im only). We adopt continuous star 
formation populations over the three timescales of 100 Myr, 1 Gyr and 10 Gyr and derive 
C70 for each timescale. The applied attenuation value is the average of the two values from 
Table 2. We then place the three model C70 values on Figure 9 (star symbols). The physical 
scales are assigned to the models based on an assumed region crossing speed of stars of v ~3 
km-s~^ (in agreement with typical velocity dispersions of young stars, Sotnikova & Rodionov 
2003), a proxy that relates physical scales to timescales. From the crude model comparison, 
one can tell that the 200 pc calibration (filled circle) is close to the 100 Myr model while 
the whole galaxy calibration (filled square) is close to a 10 Gyr model. In fact, a longer star 
formation timescale accumulates more stars with long lifetimes, which will heat the dust. 
The contribution to the IR emission from these stars will need to be removed from the SFR 
accounting, resulting in a smaller calibration coefficient for whole galaxies than for small 
regions. 
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At high surface brightness values, the spatially resolved PACS measurements tend to 
be higher than the MIPS measurements at the same wavelength, which has been attributed 
to the onset of non-linear behavior in the MIPS detectors (see, e.g., the photometric exper- 
iments by the PACS team^^'^^). If our MIPS 70 /xm measurements in the 700 pc regions 
are, indeed, slightly underestimated, we would need to increase the measured fluxes at this 
spatial scale, in order to bring them to a consistent flux scale with the 200 pc measurements. 
As a result, in order to recover the correct SFR at 700 pc, the calibration constant would 
need to be decreased, which would only increase our already found discrepancy between the 
200 pc calibration and the 700 pc calibration. In reality, we expect the effect of the discrep- 
ancy between PACS and MIPS responses at high surface brightnesses to be small on 700 pc 
scales, since the whole galaxy PACS 70 /im measurement is not significantly different from 
the MIPS one (Dale et al. 2012). Thus our estimated diffuse 70 /im fraction of the total 
emission, as derived by comparing the 200 pc scale to the integrated calibrations, is likely 
to be an accurate estimate. 

The open circle on Figure 9 is the Calibration Coefficient derived from the same com- 
bined datasets of NGC5055 and NGC6946, without the local background subtracted from 
the 70 /xm photometry. The open diamond is the Calibration Coefficient for ~700 pc re- 
gions also without the local background subtracted. Simply by comparing the filled circle 
and diamond with the open circle and diamond, we obtain about 31% and 22% of diffuse 
emission in the 70 /im IR bands when comparing whole galaxies with 200 pc and 700 pc 
regions, respectively. Thus, even when the local background is not subtracted from the 200 
pc and 700 pc apertures, we still find significant diffuse emission in the whole galaxy, and its 
fraction still increases with decreasing aperture size (Fig. 9). The persistency of the trend 
with aperture size even in the absence of local background subtraction suggests that the dif- 
fuse emission contribution in galaxies is significant in the IR bands and supports our choice 
to remove the local background in our analysis. Our choice of local background subtraction, 
by using annuli as close as possible to the aperture, is generally more robust than other 
methods (e.g., large regions backgrounds, etc.) against background variations, both real and 
spurious. Thus, we consider our results representative of the trends we should recover in 
other galaxies, as well. However, the local background subtraction is still limited by the 
resolution and the diffuse contribution associated with scales smaller than the resolution can 
not be or at least can not be effectively removed, thus producing the differences between the 
200 and 700 pc analysis. 



3ihttps://nhscsci.ipac. caltech.edu/pacs/docs/Photomctcr/PICC-NHSC-TR-034.pdf 
■^^http://hcrscheLesac.esa.int/twiki/pub/Pubhc/PacsCahbrationWeb/ExtSrcPhotoin.pdf 
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To summarize, with the high resolution Herschel PACS 70 /im images, we find an even 
larger Cjq for a smaller aperture size (~200 pc) than previous findings by Li et al. (2010) and 
Calzetti et al. (2010) with ~700 pc regions and whole galaxies, respectively, which reveals a 
diffuse 70 /xm emission fraction of about 50% beyond our 200 pc regions. We hypothesize this 
fraction is due to heating of dust by older stars not associated with current star formation. 
With the comparison to a simple model, we find this increasing trend of C70 with decreasing 
region sizes to be physically related to the associated star formation timescale, i.e. smaller 
regions related to a shorter star formation timescale have less fractional contribution from 
dust heated by stellar populations that are not related to current star formation. 



4.4. Analysis with 100 and 160 /im 

To perform the analysis at 100 and 160 /im, we also process the PACS images of these 
two wavelengths according to the methodology described in section 3.3 (Ha and Br7 images 
as well), in order to degrade all images to the same resolution of 6.9" for 100 /im and 11.5" 
for 160 /im. The adopted aperture radii are 8" for the images at the 100 /im resolution and 
12" for the images at the 160 /xm resolution, corresponding to physical sizes of about 330 
pc and 400 pc. With these aperture sizes, 65 and 47 regions for NGC 5055 and 86 and 62 
regions for NGC 6946 are selected at the 100 and 160 /xm resolutions respectively (Fig. 10). 
The final median uncertainties (and uncertainty ranges) are about 11% (10-23%) at 100 /im 
and 12% (10-26%) at 160 /im, while the Br7 measurements have similar uncertainties as 
those mentioned in section 3.3.2 at 70 /im resolution. With the extinction corrected Bv-f 
at matched resolution as the reference SFR, the correlations between both 100 and 160 /im 
and SFR are shown in Figure 11. The fit results and derived Calibration Coefficients (where 
applicable) are listed in Table 6. 

The data for NGC 6946 show a correlation with slope close to unity between the longer 
wavelengths of FIR bands and the SFR, and the corresponding Cioo and Ciqo values are de- 
rived. However, the best-fit slopes for the NGC 5055 data deviate significantly from unity, 
being smaller than 1. Moreover, the slope for the 160 /im resolution data fiattens dramati- 
cally. To understand the physical reason behind this fiattening, we use the 70 /xm emission 
as a gauge to study the behavior of longer wavelengths at different IR brightness (Fig. 12). 
For this analysis, we match the 70 /im image of the two galaxies to the 100 /xm resolution and 
to the 160 /im resolution respectively by performing the same convolution process described 
above, where each image is convolved with the PSF of the other image. Thus the resolutions 
of the matched images are about 8.7" at 100 /im and 12.7" at 160 /im, and the adopted 
aperture radii are 10" and 14" respectively. The inner and outer annuli radii are again two 
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times and four times the aperture radii. The ratios of S(100)/S(70) and S(160)/S(70) are 
plotted as a function of S(70) for both galaxies, and a linear fit through the data in log- log 
space is drawn on each panel. For the E(160)/S(70) ratio, we use the model mentioned in 
Section 4.3 (adding the 160 ^m. output) to produce a model prediction, shown as dashed 
lines on the lower panels. The difference of S(160)/S(70) ratios between NGC5055 and 
NGC 6946 is obvious, and the model prediction agrees with the NGC 6946 data. Since any 
resolution complication is taken out by matching the resolution among images, the steepen- 
ing trend of ratios for NGC 5055 is most likely coming from the fact that this galaxy is more 
inclined than NGC 6946 which increases line-of-sight overlap and confusion and/or possible 
over-subtraction of the star forming regions contributing to the local 70 micron emission. 
This is further supported by the exacerbation of the effect for increasing wavelength (lower 
resolution), which increases the line-of-sight confusion. In addition, NGC 5055 has higher 
E(100)/E(70) and S(160)/S(70) ratios on average than NGC 6946, which agrees with the re- 
sult that NGC 5055 has globally colder dust as indicated by Draine et al. (2007) and Dale et 
al. (2012), and as discussed in section 4.2. Thus, the reason for a flat trend of S(100)/E(70) 
and S(160)/S(70) ratios versus S(70) for NGC 5055, which leads to a flat trend of the 100 
and 160 /im correlation with the SFR, is summarized as: 1) a combination of decreasing 
resolution and high inclination, which produces higher line-of-sight overlap, and results in 
a higher degree of uncertainty when removing the diffuse contribution of NGC 5055 using 
local background subtraction, and/or causes over-subtraction of the star forming regions con- 
tributing to the local 70 micron emission, and/or 2) the regions in NGC 5055 have colder dust 
than that in NGC 6946, possibly a consequence of self-shielding by dust in the Hll regions 
of NGC 5055, which have high extinction, combined with a globally lower luminosity for the 
interstellar radiation field (section 4.2), and contribute more to the longer wavelengths. 

Although we derive Cioo and Ciqq from the NGC 6946 data, which has reasonably linear 
behavior, the difference between NGC 5055 and NGC 6946 questions the applicability of 
our Calibration Coefficient to general cases for 100 and 160 /xm. Conservatively, one might 
assume that the longer wavelengths are not as good SFR indicators as the shorter wavelength 
70 /xm emission, but we will revisit this issue with a larger sample of galaxies with Br7 
imaging. 



5. Summary 

We use archival CFHT WIRCam Br7 narrowband and Ks broadband observations for 
two galaxies, NGC 5055 and NGC 6946, to produce Br7 emission line maps. We adopt 
the Br7 emission as our reference SFR indicator for calibrating SFR(70) with the Herschel 
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PACS 70 /ill! data from the KINGFISH observations. The Hll region luminosity functions 
of the two galaxies have a slope consistent with -1, which agrees with previous studies using 
Ha, and lends credence to the fact that the Br7 map processing, continuum subtraction 
and absolute photometric calibration are reliable. A comparison with the free-free emission 
shows the same consistency in the Br7 maps. 

The SINGS and LVL Ha images are used for extinction correction, assuming Case B 
recombination. Both the Br7 and Ha images are convolved with the observed Herschel PSF 
of the 70 /xm band and the 70 /xm band is also convolved with the Br7 PSF, to better 
match the resolution. Sources are selected based on presence of both IR and Br7 emission 
peaks; thus only the Hll regions that are bright enough (to have detected Br7 emission) with 
enough dust (to have IR emission) are selected for our calibrations of SFR(70) and analysis. 
Aperture photometry is performed on the sources using the standard approach of employing 
annuli around the photometry apertures to remove the local background. 

We use all sources with S/N>3 at a depth of L(Br7)~10^^'^ which corresponds to star 
clusters with mass ~ 3-6x10^ Mq, assuming an average age of 4-5 Myr and a Kroupa 
(2001) stellar IMF. With the estimates of Cerviho et al. (2002) adapted to our IMF, it gives 
a fluctuation in the ionizing photon flux less than 30%. Thus random fluctuation in the 
ionizing photon flux is not a concern. Our analysis flnds: 

1. The 70 /xm emission can be combined with the Ha emission to produce an unbiased 
SFR indicator that combines dust obscured and unobscured star formation. The two 
luminosity surface densities are combined as T,Ha,obs + (0.011 ± 0.001)S(70), where 
the proportionality constant between the observed Ha and 70 /xm emission is derived 
for NGC6946 (Table 4). The same constant for NGC5055 would be higher, but this 
galaxy suffers from strong inclination effects (Section 4). 

2. Correlations between the S(70) and S(SFR) give the Calibration Coefficient estimate, 
C7o,2oopc = 1.18(±0.02) (Table 5), after combining data from both galaxies. Comparison 
between the derived C70 in this work and those from Calzetti et al. (2010) and Li et al. 
(2010) reveals a trend of increasing C70 with decreasing region size (Fig. 9), implying 
that in galaxies about 50% of the 70 /im emission outside of a 200 pc scale centered 
on bright Hll regions is unrelated to current star formation. The trend is still present 
even when the local background is not removed from the 70 /im photometry. Further 
analysis with a simple model comparison links this trend to a relation with different 
star formation timescales, i.e. larger regions have longer star formation timescales and 
thus have more IR emission produced by dust heated by stellar populations that are 
not related to current star formation as traced by the ionizing photons. This result. 
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first discussed in Li et al. (2010), is confirmed here by exploiting a larger range of 
physical sizes. 

Whenever applying a SFR(IR) calibration, it is essential to choose a suitable SFR(IR) 
at the relevant physical scale (star formation timescale) of the system. As shown by Li et 
al. (2010) (Equation (5) in the paper), a metallicity change from about solar to extreme 
sub-solar could mean a factor 2 change in the Calibration Coefficient, which should also be 
taken into consideration when dealing with low metallicity systems. 

A similar analysis for PACS 100 and 160 /xm is presented in Section 4.4. Although 
the results are limited due to the smaller number statistics with lower resolution and the 
dispersions in the data are similar to or slightly larger than that of the 70 /xm, the increased 
variations between the two galaxies at longer wavelengths hint that they are possibly not 
as good SFR indicators as 70 ^m. due to more significant contributions from dust heated 
by stellar populations that are not related to current star formation and its variations in 
different galaxies. This is in agreement with the conclusions of Calzetti et al. (2010) for 
whole galaxies and Lawton et al. (2010) for Hll regions in the LMC and SMC. Ultimately, 
given the differences we find in the calibration constants of the two galaxies analyzed in 
this study, we will need to investigate a larger sample, in order to confirm the mean values 
derived in this work. We expect to perform the same analysis using a 10-fold larger sample 
in the near future. 
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Table 1. Multi- wavelength Data 



Data 


Instrument 


Central Wavelength 


PSF 


(1) 


(2) 


(3) 


(4) 


Br7 


CFHT WIRCam 


2.166 /xm 


-1.1" 


Ks 


CFHT WIRCam 


2.146 /xm 


-1" 


Ha 


KPNO 2.1m/2.3m Bok 


6573A 


~1" 


70 /iin 


Herschel PACS 


70 fxm 


5.5" 


100 /iin 


Herschd PACS 


100 /im 


6.9" 


160 /iin 


Herschel PACS 


160 /xm 


11.5" 


Free-free emission 


GET Ka-band 


33 GHz 


25" 



Note. — col 1. data name; col 2. telescope, instrument and filter; col 3. 
filter's central wavelength; col 4. FWHM of the PSF. 

Table 2. Sample 



Galaxy 


Type 


Distance 


Diameter 


i 


12+log(0/H) 






(Mpc) 


('x') 


(deg) 


(PT) (KK) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) (7) 



[NII]/Hq E(B-V) 

Foreground Regions 
(8) (9) (10) 



NGC 5055 SAbc 7.94 

NGC 6946 SABcd 6.8 



12.6x7.2 56 8.40 9.14 0.486±0.019 0.018 0.65 

11.5x9.8 32 8.40 9.05 0.448±0.087 0.303 0.44 



Note. — col 1. galaxy name; col 2. galaxy type; col 3. galaxy distance from Kennicutt et al. (2011); 
col 4. galaxy sizes from Kennicutt et al. (2011); col 5. inclination from Moustakas et al. (2010); col 6-7. 
characteristic oxygen abundances of the galaxies; the two columns, (FT) and (KK), are the two oxygen 
abundances listed in Table 9 of Moustakas et al. (2010): the PT value, in the left-hand-side column, is 
from the empirical calibration of Pilyugin & Thuan (2005); the KK value, in the right-hand-side column, 
is from the theoretical calibration of Kobulnicky & Kewley (2004); col 8. [Nil] /Ha ratio from Kennicutt 
et al. (2009); col 9. E(B-V) used for foreground Galactic extinction correction from NED (Schlegel et al. 
1998; Schlafly & Finkbeiner 2011); col 10. median E(B-V) derived from Ha and Br7 luminosity for selected 
Hii regions at 70 fiui resolution, after removal of the Galactic foreground extinction (Schlegel et al. 1998; 
O'Donneh 1994) . 
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Table 3. Archival CFHT Images 







Ks 








Br7 




Note 


Galaxy 


# 


Exposure Time 


(s) 


# 


Exp 


osure Time 


(s) 




(1) 


(2) 


(3) 




(4) 




(5) 




(6) 


NGC5055 


112 


20 




99* 




100 




RUNID:07AD86 


NGC6946 


61 


10 




61 




110 




RUNID:07BD91 



Note. — col 1. galaxy name; col 2-3. number of images and exposure time with 
the Ks broadband filter; col 4-5. number of images and exposure time with the Br7 
narrowband filter; col 6. observation information. *12 of the narrow-band exposures 
for NGC5055 have exposure times of 200 seconds, and we rescale them to 100-seconds 
exposure maps. 



Table 4. Fitting results: log[S(70)/S(Ha, obs)] vs. log[10°-^-^("") - 1] 



Fitting Parameters 



(1) 



ao 


bo 


ai 


a. 


(J 


(2) 


(3) 


(4) 


(5) 


(6) 



NGC5055 1.787(0.034) 0.804(0.049) 1.693(0.021) ... 0.23 

NGC6946 1.881(0.027) 1.026(0.077) 1.965(0.024) 0.011(0.001) 0.24 



Note. — col.l galaxy; col 2-4. fit parameters as in y=ao+box and y=ai-|-x; 
col 5. derived a assuming SHa,corr = Sna.obs + ctSyo with a=l/10^^ (see section 
4.1); col 6. 1-cr dispersion of the data around the linear fit. 
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Table 5. Fitting results: log[S(70)] vs. log[S(SFR)] 

Fitting Parameters 
ao bo ai C70 a 

(1) m (^ (4) (5) (6) 

NGC5055 42.967(0.037) 1.176(0.047) 42.829(0.010) 1.48(0.04) 0.24 

NGC6946 43.038(0.020) 1.097(0.033) 43.001(0.009) 1.00(0.02) 0.19 

Combined 43.065(0.019) 1.232(0.027) 42.928(0.007) 1.18(0.02) 0.23 



Note. — col.l dataset; col 2-4. fit parameters as in y=ao+box and y=ai+x; 
col 5. derived Calibration Coefficient, C70 = 10"^^^"^; col 6. 1-cr dispersion of 
the data around the linear ffi. 
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Table 6. Fitting results: log[S(100)] & log[S(160)] vs. log[S(SFR)] 



Fitting Parameters 



ao bo ai C a 

(1) ^^ {^ (4) (5) (6) 

100 ;um 

NGC5055 42.818(0.026) 0.771(0.028) 43.013(0.011) ... 0.20 

NGC6946 43.079(0.027) 1.091(0.037) 43.041(0.010) 0.91(0.02) 0.21 

Combined 42.987(0.018) 0.963(0.023) 43.029(0.007) 0.94(0.02) 0.22 

160 /xm 

NGC5055 42.508(0.034) 0.503(0.035) 42.938(0.017) ... 0.15 

NGC6946 42.863(0.050) 1.117(0.068) 42.793(0.017) 1.61(0.06) 0.22 

Combined 42.628(0.026) 0.688(0.030) 42.859(0.012) ... 0.20 



Note. — col.l dataset; col 2-4. fit parameters as in y=ao+box and y=ai+x; 
col 5. derived Calibration Coefficient, C = 10^^~"^; col 6. 1-cr dispersion of 
the data around the linear fit. 
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Fig. 1. — From top left to bottom right: Br7, Bvy at 70 /im resolution, 70 /im and Ha 
images for NGC6946. The green circles are the apertures we adopted for our analysis at 
70 nm. resolution. The boxes are the local regions, with relatively even background, selected 
for the comparison between the two methods of subtracting the local background (Section 
3.3.1). North is up and east is left. The image size is about 450"x450" (about 14.8x14.8 
kpc) for each panel. The uneven background of the Br7 images is greatly enhanced after 
convolving to the 70 /xm resolution. The units for the grayscale bars are, 10~^° ergs-s~^-cm~^ 
per pixel for Br7 and Ha with pixel sizes of 0.304" and for Br7 at 70 /xm resolution with a 
pixel size of 1.4", and Jansky (Jy) per pixel for 70 yuni with a pixel size of 1.4". 
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Fig. 2. — From top left to bottom right: Br7, Br7 at 70 /im resolution, 70 /im and Ha images 
for NGC5055. The symbols and units for the grayscale bars are the same as described in 
Fig. 1. North is up and east is left. The image size is about 410"x410" (about 15.8x15.8 
kpc) for each panel. The uneven background of the Br7 images is greatly enhanced after 
convolving to the 70 /xm resolution. 
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Fig. 3. — Observed Br7 emission luminosity functions. The vertical dotted lines are 5-a 
detection limit, logL(Br7)~35.7. The dashed lines are power law fits, dN/dlogL(Br7) oc 
L(Br7)", for the luminosity functions from the 5-0" detection limit to the continuous bright 
end, which are all consistent with dN/dL(Br7)ocL(Br7)~^ or dN/dlogL(Br7)ocL(Br7)~^. 
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Fig. 4. — SFR(Br7) as a function of SFR derived from thermal free-free emission at 33GHz 
(Murphy et ah 2011a). The triangles, circles and squares are values derived using electron 
temperatures of 5000, 10000 and 20000 K, respectively. The open diamonds are values 
calculated using Br7 emission without internal extinction correction, assuming a typical 
electron temperature of lOOOOK. The filled symbols show the preferred temperature for each 
region as derived in (Murphy et al. 2011a). The one open diamond in the middle of the plot 
without a corresponding filled symbol is for Enuc 4, which we have no coverage with our Ha 
map thus no extinction corrected value. The typical error bar for SFR(Br7) is shown on the 
upper left corner. The unit is Mq • yr~^ for SFR. 
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Fig. 5. — S(HQ;,obs)/S(Br7, obs) as a function of S(70). The two panels are data for 
NGC 5055 and NGC 6946 respectively, after removal of foreground Galactic extinction. The 
dashed line indicates the extinction free value of S(Hq;, obs)/S(Br7, obs) (~ 103). Error 
bars are shown for each data point. The unit is ergs-s~^ -kpc"^ for for the luminosity surface 
density. 
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Fig. 6. — Ratio of S(70) to S(HQ;,obs) as a function of 10'^-^"^(^")-l (representing the atten- 
uated fraction of Ha emission, see section 4.1). The two plots are data for NGC5055 and 
NGC 6946 respectively. Dotted lines are linear fits in log-log space and solid lines are linear 
fit with unity slope in log-log space. 
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Fig. 7. — S(H«, obs) + 0.011S(70) as a function of extinction corrected S(Hq;). The open 
diamonds and open squares are data for NGC 5055 and NGC 6946 respectively. The dotted 
hne is a hnear fits in log-log space and the solid line is the 1-to-l line. The linear fit slope 
is 1.09(0.02). The dispersion in the data around the linear fit is 0.22 dex. The unit is 
ergs • s"^ • kpc^^ for luminosity surface density. 
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Fig. 8. — S(70) as a function of S(SFR), from extinction corrected S(Br7). The three 
columns are data for NGC5055, NGC6946 and both combined respectively. In the third 
column, open diamonds are for NGC 5055 and open squares are data for NGC 6946. Dotted 
lines are linear fits in log-log space and solid lines are linear fit with unity slope in log-log 
space. The units are M0-yr^^-kpc^^ for star formation rate surface density and ergs-s^^-kpc~^ 
for luminosity surface density. 
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Fig. 9. — Calibration Coefficients of the 70 /im band (C70) as a function of region physical 
size. The plus symbol is for NGC5055 and the cross symbol is for NGC6946, using the 
values derived in the present work. The filled circle is for both galaxies combined (with 
a combined sample of 160 regions). The filled diamond is from Li et al. (2010, using 556 
regions). The filled square is from Calzetti et al. (2010, using 189 star-forming and starburst 
galaxies). For all data with filled symbols, 3-cr error bars are shown (the filled square is 
larger than its error bar). The open circle is for the two galaxies in our work combined, 
but without the local background removed from the 70 /im photometry. Similarly the open 
diamond is for the 700 pc regions of Li et al. (2010) without the local background removed. 
The three star symbols are C70 determined for continuous star formation populations with 
different timescales, using a simple model; they are placed on the plot after associating 
physical scales to their duration of star formation (see text). 
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Fig. 10. — From top left to bottom right: 100 /im and 160 /xm images for NGC5055 and 
100 /im and 160 ^m. images for NGC6946. The green circles are the apertures we adopted 
for 100 /xm and 160 fira resolution analysis, respectively. North is up and east is left. The 
image size is about 410"x410" (about 15.8x15.8 kpc) for each panel of NGC5055 and about 
450"x450" (about 14.8x14.8 kpc) for each panel of NGC6946. The convolved Br7 images 
are not shown but the uneven background is worse with lower resolution. The units for the 
grayscale bars are all Jy per pixel, while the 100 /im images have pixel sizes of 1.7" and the 
160 yum images have pixel sizes of 2.85". 
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Fig. 11. — S(IOO) and S(160) as a function of S(SFR), from extinction corrected E(Br7). 
The vertical panels are data for NGC5055 (left), NGC6946 (middle) and both combined 
(right), respectively. In the right panel, open diamonds are for NGC5055 and open squares 
are for NGC6946. Dotted lines are linear fits in log-log space and solid lines are linear fit 
with unity slope in log-log space. We notice the flattening trend of NGC 5055 with longer 
wavelength, which is mainly due to a combined effect of decreasing resolution and higher 
inclination for this galaxy than NGC 6946. The units are Mq • yr^^ • kpc^^ for star formation 
rate surface density and ergs • s^^ • kpc^^ for luminosity surface density separately. 
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Fig. 12. — The ratios of E(IOO) and E(160) to S(70) as a function of S(70). The sohd hnes 
are hnear fits through the data in log-log space. Error bars are shown for each data point. 
The ratios for NGC5055 are on average higher than the ratios for NGC6946. The dashed 
lines on the lower panels are the same expectation from the Draine & Li (2007) model, 
which agrees with the NGC6946 data. The NGC5055 ratio of 160 /im to 70 /xm is much 
steeper than the model prediction, and is most likely related to the higher inclination of 
NGC5055, since the resolutions are matched between 70 /xm and 160 /im for this analysis. 
The corresponding U, the average stellar radiation field strength defined by Draine & Li 
(2007) and adopted in the model, are about 0.94 to 3.07 and 1.50 to 4.97 for the range of 
S(160)/S(70) ratio for NGC5055 and NGC6946 respectively. The unit is ergs ■ s"^ • kpc'^ 
for luminosity surface density. 



